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Abstract 

We consider a system of two interpenetrating Bose-Einstein condensates of 
atoms in two different hyperfine spin states. We show that in the presence of 
a small coupling drive between the two spin levels, there exist domain walls 
across which the relative phase of the two condensates changes by 2it. We 
give the physical interpretation of such walls. We show that the wall tension 
determines the force between certain pairs of vortices at large distances. We 
also show that the probability of the spontaneous decay of the domain wall is 
exponentially suppressed, both at finite and at zero temperature, and deter- 
mine the exponents in the regime of small Rabi frequency. We briefly discuss 
how such domain walls could be created in future experiments. 

I. INTRODUCTION 

One of the long-sought goals in low-temperature physics is the creation of two inter- 
penetrating superfluids. Early efforts were directed at mixtures of helium isotopes. More 
recently, following the experiments with Bose-Einstein condensates (BEG) of atomic gases 
considerable efforts have been made to create systems where two species of atoms con- 
dense simultaneously. This goal was partially achieved for two different hyperfine spin states 
of ^^Rb, which were condensed in the same trap by the technique of sympathetic cooling [Q. 
Later the dynamics of the complex relative motion of the condensates has been studied 0. 
The possibility of the measurement of the relative phase between the two condensates has 
also been demonstrated 0. In these experiments the two condensates have a substantial 
overlap, although they do not completely interpenetrate each other in the stationary state. 
A similar state called "spinor condensate" has been observed for sodium gas [0]. 

Theoretical investigation of two-component Bose systems has started many decades ago 
with the construction of the phenomenological hydrodynamic equations in the spirit of the 
Landau-Khalatnikov two-fluid model for the one-component BEG Later, this construc- 
tion has been put onto a microscopic basis p|-p!T|. Recent experiments with alkali atoms 



have revived the interest in the theory of such systems. Hartree-Fock theory has been suc- 
cessfully tested on the two-component '^^Rb system [|I^]. The stability |T^, ground-state 
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properties |]T4| and collective excitations have been studied theoretically by using the 
Gross- P it aevskii equations. 

Many properties of two- component, or binary, BEG can be understood from symme- 
try arguments. Gompared to one-component Bose superfluids, two-component systems have 
more interesting pattern of symmetry and symmetry breaking. Gondensation in binary Bose 
systems corresponds to the spontaneous breaking of two (instead of one) global U(l) sym- 
metries. These symmetries are related, by Noether's theorem, to the separate conservation 
of the number of atoms of each of the two species. The quantum state of the binary Bose 
system, therefore, is characterized by two phases of the two condensates. Gorrespondingly, 
the physics of binary BEG is also richer than of usual one-component systems. 

The effects of a symmetry are often best exposed by violating the symmetry explicitly 
in a controlled fashion. A very interesting feature, specific to systems consisting of atoms 
of the same isotope in different spin states, is that it is possible to couple two condensates 
by a driving electromagnetic field tuned to the transition frequency. In this case atoms can 
be interconverted between the two spin states and the numbers of atoms of each species 
are not conserved separately anymore; only the total number of atoms is constant. This 
implies that, in the presence of the coupling drive, only one U(l) symmetry remains exact, 
the other one is explicitly violated. The preserved U(l) symmetry obviously comes from 
the conservation of the total number of atoms, and corresponds to changing the phases of 
the two condensates by the same amount (i.e., leaving the relative phase unchanged). The 
violated U(l) corresponds to changing the relative phase between the two condensate. The 
presence of the coupling drive lifts the degeneracy of the ground state with respect to the 
relative phase. 

In this work, we show that a sufficiently small violation of the U(l) symmetry corre- 
sponding to the relative phase leads to the existence of a nontrivial static configuration — 
a domain wall inside which the relative phase changes by 27r. This configuration is a local 
minimum of the energy. However, the domain wall is not topologically stable and can "un- 
wind" itself. To unwind, however, the system must overcome an energy barrier. Thanks to 
this fact, the rate of the spontaneous decay of the domain wall is exponentially suppressed. 

Our paper is organized as follows. Section |I| introduces the field-theoretical description 
of binary BEG. In Sec. |T| we describe the domain wall configuration, whose physical inter- 
pretation is given in Sec. |IV|. Section |V| deals with the boundary of finite domain walls and 
the related phenomenon of "vortex confinement" . Section ^ contains concluding remarks. 
In Appendix ^ we find the metastability condition for the domain wall in the particular 
case when the densities of the two components are equal, and in Appendix |B| two different 
mechanisms for the decay of the domain wall, operating at different temperature regimes, 
are considered. 



II. THE LAGRANGIAN, ITS SYMMETRIES AND NORMAL MODES 

In this Section, we use field theory to describe general properties of binary BEG. Our 
goal is to introduce notations and the formalism to lay the ground for the discussion of the 
domain walls in the next Section. 

A binary dilute Bose system is described by a quantum field theory of two complex scalar 
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fields V'l and These fields have the meaning of the wave functions of the two condensates. 
The dynamics of these fields is governed by the following Lagrangian, 

L = th{4dtA + 4dt^2) - H{^,, ^2) , (2.1) 

where the Hamiltonian H{ipi,ilj2) has the form 

2m 

+ + 9i2\H'\H' - lm4^2 + 4^1) . (2.2) 

In Eq. ( p. ID /ii^2 are the chemical potentials of the two species,Q gij is the scattering ampli- 
tude, in the zero momentum limit, between an atom of the i-th species and that of the j-th 
species, and are proportional to the scattering lengths aij, 

gij = inh'^aij/m; (2.3) 

and Q is the Rabi frequency arising from the coupling drive. 

By varying the action S = JdtdxL with respect to ipi^2, the familiar Gross-Pitaevskii 
equations are directly obtained: 

h'^ 1 
ihdtiJi = ~ A^i^i + (5'iilV^iP + 9i2\i'2\'^)i'i - -^^i'2 ; 

1 

ihdtij2 = -^^^2 - f^2i^2 + igi2\H^ + g22\H^)i^2 - ^hn^Ji . (2.4) 

Let us start by finding the ground state when the coupling drive is off, f2 = 0. In the 
superfluid ground state, both ipi and ip2 have nonzero expectation values. These can be 
found by minimizing the potential energy part in Eq. (|2.1|) with respect to ipi and 1^2- This 
minimization procedure gives the equations determining the densities ni = {ipil"^, n2 = |^2p 
in terms of the chemical potentials fii and fj,2, 

guui + gi2n2 = /ii , 

fi'12^1 + 5'22^2 = fJ'2 ■ (2.5) 



More conveniently, one could view Eq. (|2.5| ) as the equations fixing the chemical potentials 
for given values of the densities. 

Strictly speaking, Eqs. (|2.5| ) only correspond to an extremum of the potential energy. For 
it to be a local minimum, the quadratic form gunf + (722^2 + 25'i2^i't-2 needs to be positive 
definite: 

9n922 - 5-12 > . (2.6) 



In real experiments ^1 and fi2 are functions of coordinates. We assume here that the trapping 
potentials are sufficiently wide so that these chemical potentials can be put to constants. 
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In fact, Eq. (|2.6|) is the condition for the mixture of the two Bose superfluids to be thermo- 



dynamically stable against segregation |T^,|T2[. In this paper we shall assume that ( p.6|) is 
satisfied. 

In principle, the constants gu, gi2 and (722 can be arbitrary. In this paper we limit 
ourselves to the regime when all three scattering lengths are close to each other, an ~ ~ 
022- In the case when the two species are rubidium atoms in different hyperfine states, these 
lengths were found experimentally to differ by no more than a few percent. The assumption 
that Qii ~ gi2 ~ 5^22 also introduces considerable technical simplifications in our treatment. 
We introduce the "average" scattering amplitude 

9 = ^ , (2.7) 

and the deviations from the average 

59 = -{912- 9) , Sg' = gu - g = -{g22- 9) , (2-8) 



so that 6g ~ 6g' <^ g. The stability condition (|2.6| ) implies that 6g > 0. Analogously, we 
introduce the average scattering lengths a and the deviations 6a and 6a'. Note that in the 
limit gn = gu = (722 the Hamiltonian has an SU(2) symmetry which leads to interesting 
imphcations jTBI . 

With the Lagrangian at hand, the discussion of symmetry in the Introduction can be 
made concrete. In the absence of the coupling drive, f2 = 0, the Lagrangian possesses 
a U(l) X U(l) symmetry with respect to independent phase rotations of the fields, 

7/^1 ^ e^^^V'i, V^2^e*"¥2. (2.9) 

The corresponding conservation laws are those of the numbers of particles of each species, 
Ni = Jd^xipl^pi and N2 = J(Pxip2ip2- That Ni and N2 at Q = are conserved separately is 
actually a basic assumption made when we wrote down the Lagrangian (|2.1| ). This assump- 
tion is not automatically satisfied: it requires that only elastic scattering between atoms is 
allowed; inelastic scattering is forbidden. For binary Bose systems made of rubidium atoms, 
this appears to be a good approximation [P!7| . 

Once the coupling drive is turned on {Q ^ 0), the Lagrangian is invariant only under a 
subset of the original U(l) x U(l) rotations; namely, those which rotate both tpi and tp2 by 
the same angle, 

V-i-^e^^i, ^2^e^>2. (2.10) 

Therefore one of the U(l) symmetries the system enjoyed at = is explicitly violated. 
Applying the Goldstone theorem, we conclude that, at = 0, there are two gapless excita- 
tions and only one of these modes remain gapless at Q ^ 0. The gapless modes at Q = are 
the phonons of the two types of sounds. One corresponds to the ordinary density wave (B 
mode in our paper, see below), and another to the concentration wave (A mode) in which 
the densities of the two species oscillate relative to each other in such a way that the total 
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density remains constant .0 If we view the two components as being made of the same atoms, 
but in different hyperfine levels, then the mode A can be alternatively interpreted as a spin 
polarization wave, in which the density of nuclear spin is a function of space and time. 

When the coupling drive is on, only the density wave remains gapless; the phonon of the 
concentration (spin density) wave is gapped. 

Let us compute the sound speeds at Q = 0. To this end we write ipi and iIj2 as 



V'l = + 5ni e^'^S = V^2 + e'^^ (2.11) 

and expand Eq. ( |2.1| ) to second order of 6n and Vy? (we will see that 5n ~ V(p). We find 



= -n{6nidtipi + 6n2dtip2) - — (ni|V(/?ip + rialVf/^sH 

2m 

--{gu5nl + 2gi2Sni5n2 + g225nl) (2.12) 

(we have thrown away total derivatives). The density fluctuations Srii can be "integrated 
out" and replaced by the saddle point values Srii = —hJ2j{g~^)ijdt'fj- As a result, Eq. (p.l2|) 
becomes 

^(2) ^^Y^^g-i^.^Q^^^Q^^^ _ ^^n,\V^,\' + n2\V^2\') . (2.13) 

Thus, the dispersion relations for the phonons are linear, u = uk, and the sound speed u 
satisfies the characteristic equation 

detf(9-'L-:if"' MUo. (2.14) 



^] mu 



n2 



When gii ^ gi2 ~ g22, the solutions are 



25, n^^ g_n_ ^^.15) 

m n m 

We see that when 5g <^ g the speed of the concentration wave (A) is much smaller than that 
of the density wave (B), ma <C mb- The system is "stiffer" towards density fluctuations than 
towards fluctuations of concentration. That modes A and B are indeed concentration and 
density fluctuations respectively is seen from the corresponding eigenvectors. The A mode 
corresponds to such fluctuations in which 



^In the hydrodynamic regime these two modes would correspond to the first and the third sounds 
respectively, according to the established classification Q. In the zero-temperature limit the density 
wave is termed Bogolyubov (zero) sound, and there is no estabUshed term for the concentration 
wave. 
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while the B mode corresponds to 



Therefore, in the A sound rii and n2 fluctuate in such a way that the overall density remains 
constant {6ni + 6n2 = 0), while the B sound corresponds to density waves in which ni/n2, 
or concentration, is unchanged {6ni/ni = 6n2/n2)- The Lagrangian (|2.13|) , in terms of the 
normal modes 

'fA = ^i-'f2 and = (2/n)(niv?i + ra2V?2) , (2.18) 

has the form 

^(2) ^ ^ r4n^ \ul\d,^Ar - (V^a)^] + [V(5.^b)^ - (V^b)^] j . (2.19) 



When the couphng drive is on, 7^ 0, one should add to Eq. ( ^.13[ ) the potential energy 



term TiQ.^nin2 cos(v9i — V52).Q The B sound is not affected by this term since it corresponds 
to such fluctuations in which ^pi = ip2- The phonons of the A sound acquire a gap 



/\ = ^2hn- 5gn[-j^J (2.20) 
for small values of Vl. 

In our further discussion we shall need the formulas for the healing, or correlation, 
lengths, which are defined via the response of the system to a static source coupled locally 
to the particle density, /ii(a;)|^/'ip + /i2(aj)|^2p. As in the case of the sounds, there are also 
two healing lengths, 

^ _ 1 Ti / y/2_ 1 (aV'^f y^'^ 

a/8 ^/m5gn \n1n2) y/2\SgJ \n1n2J 

Cb = = / . (2.21) 

2 y/mgn ^JVo■nan 

As seen from Eq. (|2.21| ), ^ ^b. This is because ^b is the correlation length of fluctuations 
of the overall density, while ^a is the correlation length of the relative density. 

If we take the values typical for present day experiments: n ~ 10^^ cm~^, a ~ 50 A, and 
bgjg ~ 10-2, then ^b ~ 0.2 /xm, and ^a ~ 3 /im. These lengths are smaller than the typical 
system size in experiments. 



^ In the Hamiltonian approach, the effect of the coupling on the spectrum can be treated by 
diagonalizing Hamiltonian using an exact transformation [p^ ]. 
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III. THE DOMAIN WALL OF RELATIVE PHASE 



The existence of the domain wall of relative phase in binary BEC with a coupling drive 
can be shown in a rather simple way. Let us first focus only on fluctuations of the fields on 
length scales much larger than the largest healing length ^a- In this case, the amplitudes rii 
and n2 of iIji^2 can be regarded as "frozen" and the only important degrees of freedom are 
the phases, i.e., (pi and ip2- The energy of the system is a functional of (pi and (p2, 



E[(pi,(p2] = J (fx 



— (ni(V<y5i)^ + n2(Vv22)^) - hny/nrn^ cos{lpi - LP2) 



(3.1) 



The potential energy term —hQy/nin2 cos(<y9i — (f2) has its minimum at (fx = fi — (f2 = 0. 
The configuration in which (fA = over the whole space is obviously the global minimum 
of the total energy, and hence is the ground state. 

The domain wall solution that we shall describe is, in contrast, a local minimum of the 
energy. To find the profile of the domain wall, we vary Eq. (|3.1|) with respect to (pi and ip2 
and obtain the following equations 

— niV^V^i = n2V'^'P2 = 'hVLJnin2 sva.{ipi — ip2) ■ (3.2) 

m m 

For domain walls, ip>i and ip2 are functions of only one coordinate, say, and aX z = ±00 
both approach constant values. A nontrivial solution to Eq. (|3.2| ) satisfying these conditions 
is 

n2 rii 
n n 

where 93A = 4arctane , k =— . (3.3) 

h y/n{n2 

The characteristic width of the domain wall is k^^. The tension of the domain wall is 

^ = 8— -77^ — ^ • (3.4) 



m 



1/2 \ ^2 I 



The relative phase ip^ changes from to 27r as z runs from —00 to +00. Note that ipA is 
defined modulo 2n so it goes a full circle as one passes through the wall. Therefore from the 
point of view of the energy functional (|3.1| ) the domain wall (|3.3|) is a topologically nontrivial 
configuration, which can not be continuously deformed into the ipA = configuration. In 
fact, one can prove that the domain wall is a configuration with minimal energy defined 
in Eq. ( p.l|) among those where ipA changes by 27t from z = —00 to z = +00, and hence, 
cannot decay away, as long as Eq. ( \3. 1\ ) applies. The domain wall is similar to the soliton 
of the sine-Gordon model. There is a small difference: the ground states on two sides of 
the domain wall are different at z = —(yo Lpi = (p2 = 0, while the ground state ai z = +00 
is (pi = 27rn2/n, 1^2 = —'iixni/n. Since 'p>i^2 are defined mod 27r the latter is equivalent to 
V?! = V52 = 11^2! n. 

In reality, Eq. ( p.l|) is not the full Hamiltonian of the system: it is only an effective 
description valid at length scales larger than both healing lengths. The full theory (p.2|) 
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contains, besides ipi and (p2, also the density fluctuations Sni and Sn2. As a consequence, 
in the full theory, the domain wall (|3.3|) can be continuously deformed into the trivial 
configuration (/^a = 0. Such deformations necessarily pass through field configurations where 
either ni or n2 vanishes at some points: at these points (p^ is ill-defined. Thus, the domain 
wall is not truly topological and can "unwind", i.e., decay away. The fact that the ground 
states on the two sides of the wall are different does not prevent the decay: it is possible to 
construct a field configuration interpolating between (pi = ip2 = and ipi = ip2 = 27m2/n 
with arbitrarily small energy per unit area. 

Although the domain wall is not a global minimum of the energy functional, it can still 
be a local minimum. In this case, to deform the domain wall into a "topologically" trivial 
configuration with = one has to overcome an energy barrier. The wall is in this case 
metastable. From our previous discussion one can conclude that, roughly speaking, the wall 
is metastable when Eq. ( p.l|) applies and is unstable when Eq. (|3.1| ) is not applicable. For 
Eq. ( |3.1| ) to be valid the wall has to be wider than the largest healing length, ^a- Since the 
width of the wall decreases as one increases the Rabi frequency Q, the wall is metastable 
only when Q is smaller than some critical value Qc- Let us define as the value of the 
Rabi frequency at which the width of the wall k~^, as defined in Eq. ( |3.3| ), is equal to to the 
longer healing length ^a in E<1- ( ^.21| ): 

= 86gn(^^^Y\ (3.5) 

The wall is metastable when Q is less than some critical value Qc of order Qq, 

n^r^no. (3.6) 

To find the exact value of flc one needs to perform a more refined stability analysis. We 
present such an analysis for the special case ni = n2 (i.e., when the densities of the two 
species are equal) in Appendix 0. Parametrically, the result is consistent with Eq. (p. 61), 
and the ratio flc/^o is found to be 1/3. 

Using the numerical values typical for experiments with Rb: ~ a few 10^^ cm~^, 5a ~ 1 
A, one finds fio ~ 100 Hz. Therefore to have a stable wall the Rabi frequency needs to be 
smaller than about 100 Hz. The domain wall cannot be thinner than the longer correlation 
length ^A, which was estimated above to be a few /im. 

Even when fl < Qc, a metastable wall can still spontaneously decay (burst). Such a decay, 
as we have said, requires overcoming a potential barrier. At sufficiently large temperature, 
the mechanism for the decay is thermal activation (over the barrier). At zero or small 
temperatures, the mechanism of the decay is quantum tunneling under the energy barrier. 
The decay through thermal activation, which is relevant for the temperatures achieved in 
current experiments, is considered in Appendix where it is shown that the decay rate is 
exponentially suppressed, since a global (macroscopic) fluctuation is required, unless fl is 
very close to Qc- 



IV. PHYSICAL INTERPRETATION OF THE DOMAIN WALL 



The domain wall solution found in Sec. |T| allows an interesting physical interpretation 
which suggests a possible way for their creation in experiments We first note that in a 
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max Vj 2 



max 2^1 max 1—2 



FIG. 1. The domain wall as superfluid flows. The two components flow in opposite directions 
perpendicular to the wall. The velocities are maximal at the center of the wall. The particle 
number of each type is not conserved due to interconversion which occurs inside the wall, with 
maximal rates on two sides of the center. 

BEC the superfluid velocity is proportional to the gradient of the phase. In a two-component 
BEC, there are two such velocities. 

Vi = — VVi , V2 = —\/ip2 ■ (4.1) 

m m 

Equation ( |3.2| ), in particular, implies that the total particle number current vanishes, niVi + 
n2V2 = 0. Individually, however, the particle number current of each species j\ = riiVi and 
= ^2^2 are nonzero. From Eq. ( p.3| ) we find 

n2 2/c _ ni 2/c _ 

^1 = TTT^ ' '"2 = T-jr^ ■ (4-2) 

n cosh Kz n cosh kz 

Thus the domain wall is a configuration where the two components flow in opposite direc- 
tions. The flow is illustrated in Fig. |l|. The velocities of the components are largest at the 
center of the wall {z = 0) and decrease as one moves toward the edge of the wall. The flow 
is concentrated on the wall; outside the wall {\z\ ^ k~^) there is essentially no flow. 
Equation ( |3.2[ ) can be rewritten in terms of the currents as 

V- Ji = -V- J2 = ny/nin2sm(fA- (4.3) 



For a stationary conflguration as the one we are considering, Eq. ( |4.3| ) means that the number 
of particles in each species is not conserved. It also implies that there is a conversion of 
atoms between the two energy levels due to the coupling drive. In the left half of the wall 
(fiA < and there is a conversion of atoms of the second type to atoms of the flrst type. 
In the right half <^a > and the conversion goes the opposite way (Fig. |I|). The rate of 
conversion is 

siiiti kz 

Q^/r^n^sm ipA = -2VL^nin2 . , n , . (4.4) 

cosh kz 

As is the flow, the conversion rate is also maximal near the wall. Far from the wall {\z\ 3> 
k~^) there is essentially no conversion. The conversion rate ( [4.4|) changes sign at 2; = 0. 
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Since different species correspond to different energy levels of the atom, energy is absorbed 
in one half of the wall and released in the other. 

The interpretation of the domain walls given above suggests a possible method for their 



creation in experiments |T9|. One starts with the coupling drive off (f2 = 0) and prepares a 
state where the two condensates flow in opposite directions (for example, by manipulating 
the traps). In such a state the relative phase (^a = v'l — is a linear function of the 
coordinate along the direction of motion (say, z). One then slowly increases VL. The domain 
walls will be created and the centers of the walls are located at the points where was an 
odd multiple of vr (ivr, ztSvr, ±57r, etc.) before VL was turned on. By changing the velocity 
of the initial relative motion of the condensates and the flnal values of f2 one can change 
the separation between the domain walls and their width. Such controlled creation of the 
domain walls, hopefully, can be achieved in future experiments. 



V. VORTEX AS THE BOUNDARY OF THE WALL AND VORTEX 

CONFINEMENT 

So far we have always considered inflnite domain walls which have no boundary. It is 
also interesting, and perhaps more realistic, to consider domain walls with a boundary. We 
shall show that the domain wall can be bounded by a vortex line. 

Suppose we have a flnite domain wall whose boundary is a closed contour C (Fig. 
We shall assume that the length of C is much larger than the width of the wall so one 
can view the domain wall as an inflnitely thin membrane stretched on C (we shall call this 
picture the "thin- wall approximation"). Let us now take another, smaller contour which has 
a nontrivial linking with C {T> in Fig. As one goes along "D, one crosses the membrane 
once, so the relative phase <^a changes by 2tt. This is exactly what one expects from a 
vortex. Therefore, C can be a vortex line. We recall that the size of the core of the vortex 
is the healing length ^a, which is smaller or of the same order as the width of the wall. 
Therefore, in the thin wall approximation, we have an infinitely thin membrane bounded by 
an infinitely thin vortex line. One should note that such a bounded domain wall will tend 
to shrink to reduce its energy, which comes from the wall tension and the tension of the 
boundary vortex. 

Conversely, for small non-zero fi, a vortex must have a domain wall attached to it to 
minimize the energy due to the nontrivial phase ip^ winding around it. This means that the 
energy of a single vortex per unit length is increasing linearly with the size of the system 
in the transverse direction. This is in contrast to the situation at f2 = when the vortex 
tension has only a logarithmic dependence on the size of the system. 

Note that there are two types of vortices in binary BEG. Those of the first type, which 
we shall call the (pi vortices, have the condensate ipi vanishing at the vortex center, while 
11)2 is nonzero. Analogously the vortices have ip2 = ^ and ipi ^ Q dX their centers. As 



^ Recalling that, in current experiments, atoms 1 and 2 correspond to different nuclear spin states 
of the same atom, we can view the wall as a stationary configuration with a flow of (nuclear) spin 
across z = plane accompanied by spin flips on both sides of it. 
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FIG. 2. The boundary of a finite domain wall. 

one goes around a (pi vortex, ipi changes by 2n, while does not change, and vice versa 
for a vortex. Thus (pA changes by either 27r or —271 for the two types of vortices, so the 
domain wall can be bounded by a vortex of either type. 

In contrast to an individual vortex, a pair of ipi and vortices, placed parallel to each 
other, will have energy per unit length which is only logarithmically divergent. That is 
because the ipA "charges" of the two vortices cancel each other, so ipA is trivial at spatial 
infinity (no winding). The same situation occurs for a pair of parallel vortices of the same 
type {ipi or 1^2) with opposite winding (however, such vortex-antivortex pair can annihilate, 
while a ipiip2 pair cannot). In a certain sense, one can talk about the phenomenon of 
"vortex confinement" : vortices exist only in pairs. This confinement should, in principle, be 
observable experimentally, in a rotating two-component BEC, which can be already created 
in a laboratory The vortices are usually identified by the density depletion at their 

cores, but can be also seen as dislocations of interference fringes due to phase singularities 



21| . With the coupling drive off (fi = 0) such a system contains an equal number of (pi and 
(p2 vortices which are distributed in space with no particular correlation between ipi and (p2 
vortices. As one turns on Q, the vortices will start to pair up and at some point the system 
will become a collection of composite objects, each being a bound state of a ipi and a ip2 
vortex (Fig. H). 



© @ ©@ 
® © ® @® 



Q.=0 0<Q.<Q.c 
FIG. 3. Vortex pairing (confinement). 

The phenomenon of vortex confinement is very similar to that of quark confinement 
in the theory of strong interaction (quantum Chromodynamics) . Similarly to our vortices, 
quarks and antiquarks do not exist as individual objects, but are confined into composite 
objects — hadrons. The analogy with quantum Chromodynamics actually stretches further. 
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If one places a <y9i vortex and a <y92 vortex at a distance larger than k~^, then a domain wall 
that connects these two vortices will be formed. The tension of the domain wall is the 
force, per unit length, that attracts the two vortices. The attractive force between the 
two vortices is thus independent of their separation, given that the latter is larger than 
the width of the domain wall. This is analogous to the confining force between a quark 
and an antiquark, which is also constant at large distances. A confinement model which 
resembles most the confinement of the vortices is the three-dimensional compact quantum 
electrodynamics considered in Ref. [^. In this theory the worldlines of electrically charged 



particles are analogous to the vortices in BEC. 

One can also imagine a system of two vortices which rotate around each other so that 
the confining force (from the domain wall) balances the centrifugal force. Such a system is 
analogous to the high-spin meson states in hadronic physics where a quark and an antiquark 
rotate around each other. 



VI. CONCLUSION 

We have shown that in a system of two interpenetrating BEC with a coupling drive, there 
exists a domain wall solution. The relative phase between the two condensates changes by 
2n as one travels through the wall. The wall solution is formally similar to the kink in the 
sine-Gordon field theory, yet it is not topologically stable and can decay. In this respect, 
the wall is more similar to a soap film, which can spontaneously burst. 

From the mathematical point of view, the domain walls discussed in this paper are similar 
to the ones which have been studied in particle physics. Such domain walls appear at least in 
two contexts: in the theory of the hypothetical axion [^] and in high-density quark matter 



The similarity is that the domain wall solution arises from the spontaneous breaking 
of an approximate U(l) symmetry. In all cases the domain wall exists only if the explicit 
violation of the U(l) symmetry (determined by the value of the Rabi frequency in our case) 
is small enough. The decay of the wall in all the examples occurs via hole nucleation. 

As to the experimental realization of the domain wall, one could be optimistic since the 
estimated critical value of the Rabi frequency, for densities and scattering lengths typical for 
the rubidium gas in recent experiments, is of the order 100 Hz which is not too small. The 
width of the wall, which might be as small as a few /im if the Rabi frequency is not much 
smaller than critical, can be also accommodated inside condensates of the size characteristic 
of present-day experiments. The apparent immediate obstacle is still the creation of a system 
where two BEC truly interpenetrate. In the recent experiments with trapped atomic gas the 
region of overlap between the two BEC is still small. One can hope the technical problems 
of making a genuine two-component BEC to be solved in the near future which would enable 
one to study the domain walls experimentally. 



^Another case discussed in particle-physics literature, domain walls in zero-density QCD [25|, is 
beyond theoretical control. 
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APPENDIX A: WALL'S LOCAL STABILITY, ni 



n2 



In the case when the densities of the two components are equal, rii = n2, the maxi- 
mal frequency Qc, below which the wall configuration is still locally stable, can be found 
analytically. For illustrative purposes we shall consider this particular case in details. 

Let us recall that when the scattering lengths an, ai2, and 022 are approximately equal 
(assuming 011022 — 012 > 0) there are two healing lengths and ^b- The healing length 
related to fluctuations of the overall density is much smaller than the one related to 
fluctuations of the relative density -C ^a- We shall be interested in the case when the 

width of the wall is much larger than ^b (but we shall not presume any relation between 
and ^a)- In this case, the total density n = iV'iP + |V'2p can be considered as frozen, 
and the system can be described in terms of three variables: 6, ipi and ip2, 




n 



(Al) 



sin(^/2)e*'^2 J ' 

where 6 runs from and vr. In terms of these three variables, the Lagrangian has the form 
L = -hn[cos'^{e/2) dt^i + sm^{e/2) dt^2] - H , (A2) 



where H is the following Hamiltonian 



H{9,ipi,ip2) 



l(V^f + cos^^(V^i)^ + sin^^(V^2f 



2m 

1 1 
— dgv?' sin^ 6 TiVtn sin 6 cos(v9 



V2) 



+ -{Sfin — 6g'n^) cos 9 



(A3) 



In Eq. ([A3|) Sfx = jji — jj2- The ground state is found by minimizing the potential term in 
Eq. (|A3|) . From here on we consider the special case when in the ground state the density of 
atoms of the two species are equal, ni = n2 = n/2 (or equivalently 9 = 7r/2.) This requires 
Sfi = 6g'n. In this case H possesses a discrete symmetry with respect to replacing 9 ^ n — 9, 
(fi ^ (f2- This symmetry is what makes it possible to find flc analytically. 

It is more convenient to use, instead of (pi,2, the normal modes <fA,B defined in Eqs. 
( ^.18| ), which in the case rii = n2 = n/2 have the form 



^1-^2, 



(A4) 



In terms of these variables the Hamiltonian becomes 



13 



H 



h n 



— -6gn sin 9 — -hQn sin 9 cos ip a 



(A5) 



In order to find the domain wall configuration, we need to extremize the energy with respect 
to 9 and (Pa,b- Varying with respect to ipB, one finds 



V- (Vv9b + cos 6'Vv9a) = . 



(A6) 



This equation determines for given 9 and ipA- The task of solving Eq. ( |A6|) becomes 
much simpler if one assumes that all variables depend only one coordinate z. In this case 



dz^B + COS 9dz(pA = 



(A7) 



which can be trivially solved: (fB = f^dz cos9 dz^fA- After eliminating ifB, the energy 
functional one has to minimize is 



H 



(y9f + sin^ ^(V<^a) 



-SgTi^ sin^ 9 — -hfln sin 9 cos (pA ■ 



(A8) 



It is easy to check that the following configuration is always a local extremum of Eq. ( [A8| ) : 

(A9) 

(Pa{z) = 4arctane'''', k"^ = —— . (AlO) 



Equation ([A9|) can be guessed from the symmetry of Eq. (|A8|) under 9 ^ n — 9. To see 
if the domain wall solution is a local minimum, one needs to expand H in the vicinity of 
(PI JATOl) . One writes 



9iz) = 9iz) + 9{z) , 

ipA{z) = (pAiz) + (Pa{z) . 

To the second order in 9 and ipA the Hamiltonian (|A8|) is 



6m 



25gnm 



+ r cos^A- (V^a)^ ^' 





h^n r 








8m - 



(All) 
(A12) 



One has to find the eigenmodes of Eq. ( |A12| ): if there are no negative modes then the domain 
wall is a local minimum of the energy; if there exist a negative mode then the domain wall 
is unstable. The second, 6'-independent, term in Eq. ( |A12| ) does not have negative modes (it 
has only one zero mode corresponding to the translation of the wall along the z direction) 
and does not lead to instability, and hence can be ignored. Taking into account the explicit 
solution ipA = 4arctane'^^, the first term in Eq. ( |A12| ) is 



8m 



(A13) 
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The well-known operator 



- - (AM) 
cosn kz 

has the lowest eigenvalue equal to — 4A;^, corresponding to the eigenfunction cosh~^ kz, which 
implies that H does not have a negative mode if f2 < Qc, where 



When Q > Qc, the configuration ([A10|) is not a local minimum of the energy functional: 



the domain wall does not exist. The value (|A15|) is of the same order as Qq in Eq- ( P-5D 
flc = 1^0- The numerical value | is specific for the case ni = n2] if ni 7^ n2 then f2c/^o is, 
in general, different. 



APPENDIX B: DECAY OF THE DOMAIN WALL 

As we have seen above, as long as < fic, the wall minimizes the energy of the system 
with respect to small local variations of the condensates. However, the global minimum of 
the energy is achieved when the phases of the condensates are constant in space. Since, 
as discussed above, the wall configuration can be continuously deformed into the global 
minimum energy configuration, i.e., the wall can decay. Here we shall estimate the lifetime 
of the wall due to such a decay. 

Since the wall minimizes the energy locally, such a deformation necessarily goes through 
a potential barrier. Thus the decay can only occur by a quantum tunneling through this 
barrier or, at finite temperature, by a thermal fluctuation over the barrier. In both cases, 
the decay rate is exponentially suppressed: in the first case by the WKB factor e~^^^, and 
in the second case by the Boltzmann factor 6"^=/^-^^, where Ec is the height of the barrier. 
The first formula applies at sufficiently low temperature, while the second one applies at 
higher temperatures. A crude estimate (see below) suggests that the corresponding crossover 
temperature is quite small (nanokelvins for parameters typical for present-day experiments), 
and is lower than temperatures achieved in present-day experiments. Thus we shall limit 
our discussion to the decay by thermal fiuctuation. This case is also simpler theoretically. 
The quantum mechanical decay, relevant for very low temperatures, will be left beyond the 
scope of this paper. 

Assuming the temperature is much smaller than the critical temperature at which one 
of the condensates melts (so that most atoms are still in the condensates), the rate of 
thermal activation across the barrier is F ~ e'^"^''^'^, where Ec is the height of the barrier 
at zero temperature. Since we are dealing with an infinite-dimensional configuration space, 
Ec should be understood as the energy at the lowest point of the barrier. This point is a 
saddle point (the energy has a single negative curvature direction). 

Some information about the exponent Ec/ksT can be obtained by a simple scaling 
argument without a detailed calculation. In terms of the dimensionless variables x defined 
as 
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X 



n 



X . 



(Bl) 



the energy becomes 



^3/2^gl/2 



d^xl- 



— (sin^ 9 + 2 cos 6^0 cos 9) — sin 9 cosfyji — ^92' 

4 2 ogn 



[B2) 



where 9q determines the relative condensate densities in the ground state at = 0: cos^^o = 
{rii — n2)/n. Consider the dependence on f2 at a given 9o. The energy functional Eq. (|B^ ) is 
equal to a dimensionful constant times a dimensionless functional which depends on Q only 
via the ratio fi/fio, where fio = {^gn/h) sin^^^ 9o, same as in ( |3.5|) . Thus the saddle point of 
the energy is also a function of this dimensionless ratio: 



If, moreover, we define a temperature Tq as 



(B3) 



m 



(B4) 



which is of the same order as the critical temperature of the Bose-Einstein condensation, 
then the decay rate for a given ni/n2 (or 9q) can be written as 



r ~ e 



-E^/ksT 



exp 



n 



To 



^1/2^^1/2^1/6 X 







(^^o)- 


= exp 



1 

V47r5ani/3 T V^^o 



(B5) 



The form of the function F{Q/Qo) cannot be found from scaling arguments alone. However, 
when ~ f2o one can expect F{Q/Qq) ~ 1, and then the exponent Ec/ksT is large, since 
5an^/^ < < 1, and T < Tq. 

The function F{Q/Qq) can be computed in the regime Q <^ Qq, where the saddle point 
configuration can be found in the "thin-wall" approximation. In this approximation field 
configurations are described at length scales much larger than the width of the domain wall. 
From this point of view the domain wall is an infinitely thin membrane. The saddle point 
configuration is a membrane with a round hole in it (Fig. ^). The radius of the hole R must 
be much larger than the width of the wall for the thin wall approximation to be valid. As 
discussed in Sec. |V|, the rim of the hole must be a vortex, since it is the boundary of the 
domain wall. 

There are two contributions to the energy difference Ec between the saddle point and 
the domain wall configurations. One contribution is negative and comes from the hole (since 
the hole is the absence of the wall). Another contribution is positive and comes from the 
rim. Therefore, 



E = 2'kRv - TrR^a . 



(B6) 
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FIG. 4. A hole in a domain wall. 



where v is the energy per unit length of the vortex (the vortex tension) and o is the domain 
wall tension. The energy ( |B6[ ) has a maximum when the radius of the hole is 



= v/a . (B7) 

This is the radius of the critical hole. Indeed, if a hole with a radius R > Rc is nucleated, 
then it will expand and eventually eat up the whole wall. If, in contrast, the radius of the 
nucleated hole is less than critical, then the hole will shrink and disappear. 
Substituting R = Rc in Eq. ( p6|) , we find the height of the energy barrier 

(B8) 



As there are two types of vortices, in Eq. ( P8|) z/ refers to the vortex with the smaller tension. 
The tension of a straight vortex is logarithmically divergent: 



h'^Tii „ , R 



. . 27rln-. (B9) 
2m 4a 

where the index i = 1,2 refers to the two types of vortices, and R is the long-distance cutoff. 
The role of R is played by either the size of the critical hole Rc or the width of the wall 
k~^. We shall see at the end of this Section that the two lengths differ only by a logarithm, 
which does not affect Eq. (p^). Therefore, to logarithmic accuracy, the vortex tension is 

Ui = In — - . (BIO) 

The argument of the logarithm is (/c^a)~^ = i^o/^Y^'^^ and is large when n <^ ^Iq, which 
justifies the use of the logarithmic approximation. 

Without loss of generality, we can assume that ni < Then the vortex of the first 
type has the smallest u (|B9|) . Substituting Eqs. ( |B1CI|) and ( p.4D into Eq. (p8|), one finds 
that the barrier height has the form of Eq. (Bl) , where 
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The decay rate (|B5| ) can be rewritten in the following form 



thermal 




^ In^^ , (B12) 



n 



\ 2/3 



am, 



(B13) 



— the critical temperature for the smallest of the two condensates {ni by our choice). The 
rate ( |B12| ) is exponentially suppressed when ^ fip and T ^ Td. 

To check the consistency of our assumptions, we note that the radius of the critical hole 

is 



a 16 V TTiiJ, n2\n^J is 



Comparing with the width of the wall in Eq. ( p.3| ), we find that Rc is larger than k'^ by a 
factor ln(f2o/^^), which is assumed to be parametrically large. Thus, the use of the thin wall 
approximation is justified. 

To make a good estimate of the crossover temperature, below which crossover proceeds 
via quantum tunneling, rather than via thermal activation, one needs to estimate the ac- 
tion on the tunneling trajectory and compare it to the exponent in the thermal activation 
rate Ec/kBT. For a very crude estimate we can take S ~ EcRc/ua- The exponents S/h 
and Ec/ksT become comparable at a temperature of order few nanokelvin, given typical 
parameters of present day experiments. 
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